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ABSTRACT 

New interferometric observations of SiO J = 5 — » 4 circumstellar line emission around the carbon 
star IRC+10216, using the Submillimctcr Array, arc presented. Complemented by multi-transition 
single-dish observations, including infrared observations of ro- vibrational transitions, detailed radiative 
transfer modelling suggests that the fractional abundance of SiO in the inner part of the envelope, 
between «3-8 stellar radii, is as high as wl.5 x 10~ 6 . This is more than an order of magnitude 
higher than predicted by equilibrium stellar atmosphere chemistry in a carbon-rich environment and 
indicative of the importance of non-LTE chemical processes. In addition to the compact component, 
a spatially more extended (r e ~ 2.4 x 10 16 cm) low-fractional-abundance (/q ks 1.7 x 10 -7 ) region is 
required to fit the observations. This suggests that the majority of the SiO molecules are effectively 
accreted onto dust grains in the inner wind while the remaining gas-phase molecules are eventually 
photodissociated at larger distances. Evidence of departure from a smooth wind is found in the 
observed visibilities, indicative of density variations of a factor 2 to 5 on an angular scale corresponding 
to a time scale of about 200 years. Additionally, constraints on the velocity structure of the wind are 
obtained. 

Subject headings: stars: abundances - stars: AGB and post-AGB - stars: carbon - stars: circumstellar 
matter - stars: individual (IRC+10216) - stars: mass loss 



1. INTRODUCTION 

IRC+10216 (CW Leo) is a cool red giant star located 
on the asympt otic giant branch (AGB) . It is classif ied as 
a carbon star ljMillerlll970llHerbig & Zappalall97Ch with 
a C/O-ratio of > 1.4 in its atmosphere (Winter s et al.l 
119941: lGlassgoldlfl996|) and is currently losing mass at a 
rate of « 1-2 x 10 -5 M Q yr _1 , as determined from a wide 
variety of observations from infrared to radio wavelengths 
and using molecular line as well as dust continuum emis- 
sion (e.g.. iKead v et al.l Il988t iG roenewe gen et af]ll998t 
ISchoier fc Olofssonll200lHSchoier et al.ll2002D . It is com- 
monly thought that IRC +10216 is in its last evolutionary 
stage on the AGB before the planetary nebula ejection 
phase. The high mass-loss rate of IRC +10216 makes 
direct observations of the stellar atmosphere difficult, 
and the future evolution of the star is possibly best 
followed through monitoring of its circumstellar mate- 
rial, although some stellar light emerges due to the in- 
homogeneous nature of the wind. Departures from a 
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homogeneous wind are evident from images of interstel- 
lar light scattered in the circumstellar envelope (CSE ) 
around IRC+10216 l|Ma,uron fc Huggins) IT99fl 12000). 
The high mass-loss rate of IRC +10216, coupled with its 
close proximity to the Sun («120pc), has made it the 
prime target for studies of not only carbon stars, but also 
AGB stars in general. 

Presently, there are 63 molecular species detected in 
AGB stars with about half of them detected in only 
IRC+10216 (|Olofssonll2006T) . A significant number of 
these molecules are unique to the circumstellar medium 
which makes the study of stellar winds important also for 
testing our knowledge of interstellar chemistry in general. 
Most of the abundance estimates are based on rather 
simple methods and are typically order of magnitude 
estimates. Detailed chemical modelling of IRC +10216 
does a reasonable j ob in explaining many of the ob- 
served abundances l)Millarl 12003^) . However, there are 
some notable exceptions. The dete ction of H2O toward s 
IRC +10216 by the S WA S satellite lIMelmck et al.l2001ft 
and reconfirmed by Odin observations (Hasceawa et alJ 
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Fig. 1. — Velocity channel maps of SiO J = 5^4 line emission from IRC +10216 obtained using the SMA. The contour levels are l.On 
Jy beam - 1 , where n = —4, —2, 2, 4, 8, 16, 32 (negative values have dashed contours), and the beam size is 3'.'3 X 2''0 with a position angle 
of 0?2 as indicated in the lower right corner of each panel. The velocity channels (given in the LSR frame and indicated in the upper left 
corner) have been binned to 4kms _1 . The systemic velocity is — 26kms -1 as determined from CO observations. Offsets in position are 
relative to a 2 000 = 09 h 47 m 57!39, feooo = 13°16'43'.'9. 



120061) came as a surprise. Recently, lHasegawa et al.l 
(2006) derived an abundance of NH3, based on obser- 
vations by the Odin satellite, that is four orders of mag- 
nitude higher than predicted by ste llar equilibrium chem- 
istry l|Gherchneff fc Barkerl Il99l . and processing by 
shock s does not seem to increase t he abundance any fur- 
ther l) Willacv fc Cherchnef3ll998|l . A possible explana- 
tion for these high abundances cou ld be Fischer- T ropsch 
catalytic processes as suggested bv lWillacvl {2004) or, in 
the case of H2O, eva poration of Kuiper-belt like objects 
l|Melnick et alJl20?)l . 

The first two more detailed studies of circumstel- 
lar abundances in larger samples of sources have 
been performed by iGonzalez Deleado et all l)2003j) and 
ISchnier et all ll200fl) for SiO in 45 M-type (C/O < 1 in the 
photosphere) AGB stars and 19 carbon stars (C/0> 1 
in the photosphere), respectively. Average SiO fractional 



abundances were obtained from a detailed radiative 
transfer analysis of multi-transition single-dish observa- 
tions. Interestingly, for the M-type AGB stars the de- 
rived abundances are generally much lower than expected 
from photospheric equilibrium chemistry («5xl0~ 5 , 
Duari et al. 1999). For the carbon stars, on the other 
hand, the derived abundances are on the average two or- 
ders of magnitude higher than predicted by photospheric 
equilibrium chemistry (« 5 x 10~ 8 , Millar 2003). In fact, 
when comparing the two distributions of SiO fractional 
abundances there appears to be no way of distinguish- 
ing a C-rich chemistry from that of an O-rich based on 
an estimate of the circumstellar SiO abundance alone. 
Moreover, there is a clear trend that the SiO fractional 
abundance decreases as the mass-loss rate of the star in- 
creases, as would be the case if SiO is accreted onto dust 
grains. 
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Fig. 2. — Azimuthally averaged visibility amplitudes (mean value taken over 2kms bins) obtained towards IRC +10216 using the 
SMA, as a function of distance to the phase centre in kA (also shown on the upper abscissa is the equivalent spatial resolution in "). The 
observations are compared with various models for the SiO abundance distrubution. Models 1—4 contain a single SiO abundance component 
whereas Model 5 in addition has a compact high-abundance component (see text for details). 



Further support for such a scenario comes from inter- 
ferometric observations of the tw o M-type AGB star s 
R Dor and L 2 Pup performed bv ISchoier et alJ {2004). 
In their analysis they found evidence of an inner com- 
pact component of high fractional abundance, consistent 
with predictions from stellar atmosphere chemistry. In 
addition, an extended low-abundance component, as ex- 
pected if SiO is effectively depleted onto grains in the 
inner wind, was required in order to fit the observations. 

Presented here are new interferometric observations of 
SiO v = 0, J = 5 — > 4 line emission from the circumstellar 



envelope around the carbon star IRC +10216 obtained 
by the Submillimeter Array (SMA). Supplemented by 
additional multi-transition single-dish radio and infrared 
observations the radial abundance profile of SiO is de- 
termined. In addition, information on the acceleration 
region close to the photosphere is obtained. 



2. OBSERVATIONS 
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Fig. 3. — Spectrum at the central pixel in the CLEANed map 
(histogram) overlayed by the best- fit model (Model 5; solid line). 
The velocity resolution is lkrns" 1 . 



SMA 1 observations of the SiO J = 5^4 line at 
217.10498 GHz were made on 2005 February 23 using the 
compact array configuration with sev en 6 m antennas . 
Details of the SMA are described bv IHo et all (|2004|) . 
The weather was excellent, with a 225 GHz atmospheric 
opacity of 0.05 (measured at the nearby Caltech Sub- 
millimeter Observatory) and system temperatures (DSB) 
ranging between 150 to 300 K. The projected baselines 
ranged from 11 — 50 kA resulting in a synthesized beam 
of 3'.'3 x 2'.'0, with a PA of 0?2 using uniform weighting. 
The digital correlator has a bandwidth of 2 GHz and the 
spectral resolution was 0.8125 MHz, corresponding to a 
velocity resolution of 1.1 kms . Phase and amplitude 
calibration were performed on the quasars 0851+202 and 
1058+015. Bandpass calibration was obtained on Cal- 
listo and Jupiter. Observations of Callisto provided the 
flux calibration; the uncertainty of the flux scale is esti- 
mated to be ~20%. The data were calibrated using the 
MIR software package developed originally for the Owens 
Valley Radio Observatory and adapted for the SMA. The 
calibrated visibility data were imaged and CLEANed us- 
ing MIRIAD. 

The actual analysis and comparison with models are 
carried out in the uv-plane to maximize the sensitivity 
and resolution of the data. Thus, we expect to obtain 
usable information on scales as low as 2", corresponding 
to the longest baselines. 

3. SIO BRIGHTNESS DISTRIBUTION 

The velocity channel maps of the SiO J = 5 — > 4 line 
emission obtained by the SMA, presented in Fig.^ show 
that the brightness distribution towards IRC +10216 ap- 
pears to have an overall circular symmetry and suggests 
that the emission is moderately resolved. At this level 
of resolution no signs of deviations from a homogeneous 
wind such as mass loss modulations or clumps are evi- 

1 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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Fig. 4. — SiO v = 0, J = 5 — > 4 spectra for 
IRC +10216. The solid line is the SEST single-dish observation 
IGonzal cz Dclgado ct al. 20ti3), whereas the histogram shows the 
spectrum at the phase centre derived from the SMA data convolved 
to the size of the SEST beam of 23". The velocity resolution is 
1 kms -1 . 
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Fig. 5. — Size (FWHM; panels (a) and (c)) and flux (panels 
(b) and (d)) of the observed SiO J = 5^4 line emission toward 
IRC +10216 as function of velocity channel (given in LSR frame), 
averaged in 1 kms - bins, as estimated from Gaussian fits to the 
visibilities. Panels (a) and (b) are fits to the SMA data alone while 
panels (c) and (d) have the flux fixed at the zero-spacing value. The 
observations are compared with the results from single SiO compo- 
nent models with (Model 4; dashed-dot-dot-dot line) and without 
(Models 1-3; dash-dot, dotted and dashed-line, respectively) a ra- 
dial velocity gradient and a model with an additional compact high 
abundance component (Model 5; solid line). A fit to the data using 
Eq.0for the size of the emitting region is also shown (panels (a) 
and (c); thin solid line). 



dent. 

In Fig. [21 the azimuthally-averaged visibility ampli- 
tudes, mean value taken over 2 kms -1 bins, are plot- 
ted as a function of distance to the phase centre in kA. 
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It is evident that there is some degree of patchiness in 
the distribution. The visibility amplitudes show two re- 
gions where the amplitudes flatten out: at around 20 
and 32 kA. The plateau at 20 kA is most evident at the 
blue-shifted velocities, while the plateau at 32 kA is seen 
at both blue- and red-shifted velocities. At the extreme 
velocities, these plateau regions are not apparent. The 
fact that these patterns do not appear constant from ve- 
locity channel to velocity channel make us believe that 
they represent real structure in the SiO brightness dis- 
tribution. The spatial scales involved are « 5" and rj 3" 
for the 20 and 32 kA plateau regions, respectively. From 
the radiative transfer modelling performed in Sect. 14.21 
the spread in the observed visibilities can be accounted 
for by SiO abundance modulations of a factor « 2-5, 
or a similar change of the H 2 density. Both the time 
scales, «200yr, and density modulations are consistent 
with those obtained from the s cattered light images of 
iMauron fc Hugginsl lfl999l I200H . 

The spectrum at the central pixel in the CLEANed 
map is presented in Fig. El (histogram) . The spectrum 
is characteristic of a well-resolved, with respect to the 
beam, emitting region of moderate optical depth. Op- 
tical depth effects are apparent on the blue-part of the 
spectrum in Fig. El which generally has a lower intensity 
than the red part due to self-absorption of emission in 
gas moving towards the observer. 

As interferometers lack sensitivity to large scale emis- 
sion it is of interest to ascertain the missing flux of the 
SMA data. In Fig. H the SMA SiO J = 5 -> 4 data has 
been convolved with a 23" circular beam (histogram) to 
represent the beam of the Swedish-ESO Submillimetre 
Telescope (SEST). The SEST spectrum has been con- 
verted from main beam brightness temperature scale to 
Jy using 

s = rj mh T~ 1 T mh , (1) 

where the main-beam efficiency rj m b = 0.5 and the sen- 
sitivity r~ 1 =40JyK -1 The telescope parameters are 
taken from the SEST homepage 2 . A comparison with 
the equivalent SEST spectrum from Sch oier et al"l (2006) 
(solid line) shows that much of the emission is resolved 
out by the interferometer. In addition, the CLEANed 
image of an extended source with missing short-spacing 
data results in an extended depression of negative surface 
brightness on which the source emission resides. This ar- 
tifact is most evident in the central velocity channels in 
Fig. ^ where the real envelope emission is expected to 
be the most extended. The negative feature will further 
reduce the size and flux of the emission in those chan- 
nels. In the central velocity channels of Fig.EJonly about 
10 Jy is recovered after the cleaning compared to » 29 Jy 
that the interferometer actually measures (Fig. EJ) . 

Assuming that the emission has an overall spherical 
symmetry, circularly symmetric Gaussians have been fit- 
ted to the visibilities in lkms -1 velocity bins. The re- 
sulting FWHMs and fluxes of the circular Gaussians fit- 
ted to the visibilities are plotted in Figs. Et and Et- 
There appears to be no systematic trend in the derived 
offsets with velocity. The flux in each velocity bin is sig- 
nificantly higher than the flux obtained from the full map 
(Fig. El histogram). The estimated flux from the Gaus- 
sian fit is, however, the predicted zero-spacing flux based 

2 www . Is . eso . org/lasilla/Telescopes/SEST/SEST . html 



on the observed visibilities. Hence, it tries to compensate 
for the missing flux in the observations. The fluxes re- 
ported in Fig. Et are lower than the total flux picked up 
by the SEST single-dish observation (Fig. El solid line) 
indicating a departure of the brightness distribution from 
that of a Gaussian at larger spatial scales. 

The overall variation in size with line-of-sight velocity 
(Fig. Et) is not as expected for a well resolved envelope 
expanding at a constant velocity, where gas moving or- 
thogonal to the line of sight subtends a larger solid angle 
than radially moving gas near the extreme velocities, as 
described by 

r 21 1 / 2 

R(v) = R s , (2) 

where R(v) is the estimated size at velocity v, R s the to- 
tal size of the emitting region, v\ SI the systemic velocity 
(LSR) , and v e the expansion velocity. Using Eq. El with 
i? s = 4'/7, «i sr = — 26.5kms~ 1 and v e = lSkms -1 , does 
not provide an acceptable fit (Fig. Et; thin solid line). 
The largest angular extent instead appears close to the 
blue part of the emission with 5'.'4 ± O'.'l in the velocity 
range —40 to —38 km s~ x (LSR). This corresponds to a 
radial size of 1.0 x 10 16 cm at a distance of 120 pc (see 
below). It should be noted that in this region the line 
profile is heavily affected by self-absorption as seen in 
Fig. El 

Self-absorption would naturally explain the larger spa- 
tial extent near the blue edge of the spectrum since pho- 
tons emitted by the gas close to the star, and which in 
the optically thin case would have no further interaction 
with the gas, can become absorbed and re-emitted at 
larger distances from the star thereby increasing the size 
of the brightness distribution. The lower intensity of the 
emission naturally follows from the decrease of density 
and temperature with radial distance from the star. 

However, there is some concern that the size-velocity 
relation in Fig. Et is affected by the lack of short-spacing 
information. The predicted zero-spacing flux near the 
systemic velocity from the Gaussian fits (Fig. Et) is 
about a factor of two lower than what is observed us- 
ing the SEST (Fig. H solid line). If instead the flux of 
the Gaussian fit is fixed at the SEST value (Fig. Eli), a 
different size-velocity relation is found, as illustrated in 
Fig. Et- The observed size- velocity relation is now bet- 
ter described by Eq. El in this case with R s = 6'.'0 (thin 
soild line in Fig. Et), corresponding to 1.1 x 10 16 cm~ 1 
at the distance of IRC+ 10216. It is noted that the size 
estimate near the blue edge of the spectrum is still sig- 
nificantly larger than expected from Eq. El due to the 
self-absorption. 

4. SIO EXCITATION ANALYSIS 

4.1. Radiative transfer model 

In order to determine the molecular excitation in 
the CSE a detailed non-LTE radiative transfer code, 
based on the Monte Carlo method, was used. The 
code is described in detail in Schoier & Olof ssonl 1)200 ll ^ 
and has been benchmarked, to high accuracy, against 
a wide variety of molecular-l ine radiative transfer codes 
in Ivan Zadelhoff et alJ l)2002fl and van der Tak et al. (in 
prep.). The CSE around IRC+10216 is assumed to be 
spherically symmetric, produced by a constant mass-loss 
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the excitation analysis to the adopted SiO fractional abundance 
(/o) and envelope size (r c ). Contours are drawn at the 1, 2, and 3<r 
levels. Also indicated is the lowest reduced \ 2 (x? e d = Xtot/C^ — 
2)) in the map. 



rate (M), and to expand at a constant velocity (v c ). A 
velocity gradient (acceleration region) is also tested in 
Sect.lOI 

The physical properties, such as the density, temper- 
ature, and kinematic structure prevailing in the circum- 
stellar e nvelope around IRC+10216 have been deter- 
mined in lSchoier fc Olofssonl (120001 12001^ . ISchoier et all 
(2002), and iSchftierlT'al.l $2Q0(f\ . based on radiative 
transler modelling of multi-transition CO line observa- 
tions, from millimetre to IR wavelengths. This model, 
where the circumstellar envelope is formed by a mass-loss 
rate of 1.5 x 10 _5 M Q yr _1 and expanding at a velocity 
of 14.0 km s -1 , is used as input to the SiO excitation 
analysis. 

The excitation analysis includes radiative excitation 
through the fi rst vibrational l y exci ted (v = 1) state, for 
SiO at 8 fim. ISchoier et all 1)2006(1 found that for SiO 
it is particularly important that this is tre ated correctly- 
Releva nt molecular data are summarized in lSchoier et al.l 
( 2005) and are made publicly available through the Lei- 
den Atomic and Molecular Database (LAMDA) 3 . For 
IRC+10216 thermal dust emission provides the main 
source of infrared photons which excite the v = 1 state. 
The addition of a dust component in the Monte Carlo 
schem e is straightforward as described in ISchoier et al.l 
(2002). The dust-temperature structure and dust- 
density profile for IRC+10216 are obtained from radia- 
tive transfer mod elling of the spectral energy distribu- 
tion using Dusty ( Ivezi c fe Elitzur| ^97j)^ Detail s on th e 
dust modelling can be found in ISchdier et al.l (2006). 
The best-fit model is obtained for an optical depth at 
10 /im of 0.9, a dust-condensation temperature of 1200 K 
(ro = 1.7 x 10 14 cm) and an effective stellar temperature 
of 2000 K. The luminosity is 9600 Lqj obtained from 
a period-luminosity relation l)Groenewegen fc Whitelockl 

3 http : //www . strw . leidenuniv . nl/~moldata 



1996), and the corresponding distance is 120 pc. It 
should be noted that this period-lumonosity relation was 
statistically derived for 54 carbon stars with pulsational 
periods in the range 150 — 520 days. An additional uncer- 
tainty comes from the fact that the period of IRC+10216 
is 630 days. 

4.2. The SiO abundance distribution 

The abundance distribution of SiO is initially assumed 
to be described by a Gaussian 



f(r) = fo exp - 




(3) 



where / = n(SiO)/n(H2), i.e., the ratio of the number 
density of SiO molecules to that of H2 molecules. Here, 
/o denotes the photospheric fractional abundance of SiO 
and it is assumed that some process effectively destroys 
the molecules at r > r c such as, e. g., photodissociation 
by the ambient interstellar uv-field. ISchoier et alJ ( 2006) 
found that such an abundance distribution could ex- 
plain multi-transition single-dish radio- line SiO observa- 
tions of IRC+10216 well, with f =2.8 ±1.1 x 10~ 7 and 
r c = 1.9±0.7x 10 16 cm (see Fig.®. 

Fig. |21 shows the result of applying the (u,v) sampling 
of the SMA observations to various SiO envelope mod- 
els. All models (1-3) are within the 1 a confidence level 
when compared to the constraints put by the single- 
dish data (see Fig. EJ). A model with /n=3.9 x 10~ 7 
and r e = 1.2 x 10 16 cm (Model 1; dash-dotted line) pro- 
vides a poor fit with too much flux on most baselines 
while a model with f = 1.7 x 10~ 7 and r e — 2.6 x 10 16 cm 
(Model 2; dotted line) generally has too little flux 
and gives also a poor fit. However, a model with 
fo = 2.3 x 10~ 7 and r c = 1.9 x 10 16 cm (Model 3; dashed 
line) provides a better overall fit, illustrating the useful- 
ness of the SMA data in further constraining the SiO 
abundance distribution. The models are summarized in 
Table [J 

The models are particularly poor at reproducing the 
observed visibilities in the most blue-shifted part of the 
spectrum, between — 42kms~ 1 to — 38kms . This indi- 
cates that there is too much self- absorption of the emis- 
sion in the models. This could possibly indicate that 
the medium is clumpy to some degree. Nevertheless, the 
size- velocity behaviour near the blue-shifted edge as seen 
in Fig. [5] is nicely reproduced by the model and the large 
observed size is explained by self-absorption in the gas 
moving towards the observer. 

Models with a constant expansion velocity, such as 
models 1-3, provide poor fits to the observed size- velocity 
distribution closer to the systemic velocity. The effect of 
a velocity field in explaining the observed size-velocity 
distribution and departures from the single Gaussian 
fractional abundance distribution of SiO are tested in 
Sect. 14.31 where also additional constraints provided by 
ro-vibrational transitions in the infrared are used. 

4.3. Probing the inner wind 

Infrared observations of ro-vibrational transitions have 
been shown to be a sensitive tool for probing both the 
physical propert ies of the inner wind as well as the abun- 
dance structure l|Keadv fc Ridgwavl ll993: Mo nnier et alJ 
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Fig. 7.— Calculated SiO v = l->0 J = l-»0 R(0), u = l->0 
J = 3-* 2 R(2) and u = l->0 J = 6^5 R(5) line profiles for 
the best-fit single-component model using /o = 2.3 X 10~ 7 and 
r c = 1.9 X 10 16 cm, with (Model 4; dash-dot-dot-dot line) and with- 
out (Model 3; dashed line) a gradient in the expansion velocity of 
the wind. The solid line is a two-component model (Model 5) in- 
cluding a compact (r c = 4.5 X 10 14 cm) high-fractional-abundance 
(/j = 1.5 X 10 — 6 ) SiO region in addition to the more extended 
(r c = 2.4 X 10 16 cm) low-fractional-abundance (/o = 1.7 X 10~ 7 ) re- 
gion and an accele ration region. The histog rams are the observa- 
tions performed bv lKea dv & Ridgway 119931) . 



TABLE 1 

Model summary 5 





h 


fa 




[cm] 


13 


Model 1 




3.9 x 10" 


7 


1.2 x 10 16 


0.0 


Model 2 




1.7 x 10" 


-7 


2.6 x 10 16 


0.0 


Model 3 




2.3 x 10" 


-7 


1.9 x 10 16 


0.0 


Model 4 




2.3 x 10" 


-7 


1.9 x 10 16 


0.5 


Model 5 


1.5 x icr 6 


1.7 x 10" 


-7 


2.4 x 10 16 


0.15 



a The abundance distribution of SiO is assumed to be a Gaussian 
described by fo and r e (see Eq.0. /j is an increased abundance in 
the inner region of the wind and j3 the slope of the radial velocity 
gradient (see text for details). 



2000). However, the only carbon star to have been stud- 
ied is IRC +10216. Observations of SiO ro-vib transi- 
tions (v = 1— > 0) for IRC + 10216 have been performed by 
iKeadv &: Ridewavl l)1993j) . They found that the observed 
transitions (from R(Q) up to -R(5)), exhibiting classi- 
cal P Cygni-type line profiles, can be reasonably well 
modelled with an SiO fractional abundance of 8 x 10~ 7 . 
Our best-fit single-dish model for IRC +10216, using 
fo = 2.6 x 10~ 7 and r c = 1.7 x 10 16 cm, provides too little 
flux in these observed ro-vibrational transitions as illus- 
trated in Fig. □ for the i?(0), R{2) a nd R(5) transitions 
(dashed lines; the original data from IKeadv fc Ridewavl 
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Fig. 8. — Adopted velocity field (solid line), based on a dynamical 
model, and the resulting H2 density structure (dashed line) using 
a mass-loss rate of 1.5 X 10 — 5 Mq yr — 1 . Also shown is the derived 
SiO fractional abundance distribution (dotted line). 



(1993) are represented by the histograms). In particular 
the absorption feature produced in the model appears 
too narrow and is located at a too large velocity. A way 
to remedy this and fit the data better is to introduce a 
velocity field, i.e., an acceleration region in the inner en- 
velope. Evidence for an acceleration region also comes 
from the estimated sizes from Gaussian fits in the uv- 
plane (Fig. [^J, where the models (1 — 3), which behave 
as expected for a resolved wind expanding at a constant 
velocity, give a poor fit to the size-velocity relation. The 
observations suggest that parts of the emission arises in 
the acceleration region. 

We find that by a dopting a velocity field described by 
l|Habing et alJll99^1 

,/3 



v(r) — v(oo) (l 



(4) 



with /3 = 0.5, it is possible to further improve the fit 
to the observed line shape. Eq. 0] with (3 = 0.5 is valid 
if the ratio of the dust and gas velocities is constant 
throughout the wind. When introducing a velocity gra- 
dient the density structure will also change according to 
the continuity equation, keeping the mass-loss rate con- 
stant. This has been taken into account in the mod- 
elling. All models are required to provide fits also to the 
single-dish radio-data within the la level. In Fig. [7| a 
model with r = 1.7 x 10 14 cm (adop ted to coincide wit h 
the dust condensation radius, see ISchoier et alJ {2006)) 
and ,9 = 0.5 (Model4; dash-dot-dot-dot line) still does 
not provide an acceptable fit. Increasing (5 even fur- 
ther, thereby extending the acceleration region, does not 
improve the fit, and also gives too narrow lines for the 
single-dish radio-line spectra. 

Keeping Eq. 01 with j3 as a free parameter, and fit- 
tin g it to a dynamical mo del of IRC +10216 performed 
bv iRamstedt et aD l)2006|) a value of /3«0.15 provides 
a good fit. A lower value for 8 means a smaller a ccel- 
eration region. W e note that IKeadv et al.1 l)1988l) and 
IKeadv k, Ridgwavl l)1993|) derive a more elaborate ve- 
locity field that also points to a small acceleration re- 
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Fig. 9. — Best-fit two-component model (Model 5; solid lines) overlayed on observed single-dish spectra (histograms). This model also 
provides the best fit to the SMA and ro-vibrational data. 



gion. However, their results were not based upon any 
dynamical model. It turns out that a model with a [3 
smaller than 0.15 provides a worse fit. The velocity of 
the wind is within 10% of its terminal value at distances 
>5 x 10 14 cm. 

The only way to improve the fit to the observed ro- 
vibrational transitions is to increase the SiO fractional 
abundance in the acceleration region. If a model with a 
compact (r e = 4.5 x 10 14 cm) high-fractional-abundance 
(/j = 1.5xl0~ 6 ) component is used a much better fit 
to the data is obtained (Model 5, which includes an ex- 
tended low abundance component, see below, that does 
not contribute to the ro-vibrational lines; solid line in 
Fig. 0). However, the region of high abundance can 
not be extended all the way down to the photosphere. 
The velocity field inside the dust condensation radius 
(ro = 1.7 x 10 14 cm) is assumed to be 2kms~ 1 in accor- 
dance with lKeadv et all l)1988|) . We emphasize that this 
is just a test and that the structure of the velocity field 
close to t he photospher e certa inly is more complex as 
shown by Wi nters et alJ <)2000j) . A much lower value of 
« 3.0 x 10 -8 is required in this region otherwise too much 
absorption will result in the ro-vibrational lines at veloc- 
ities <2kms~ 1 . The adopted velocity field, resulting H2 
density structure from the continuity equation, and the 
derived SiO fractional abundance structure are shown in 
Fig. El 

An extended low abundance component is still needed 
in order to model the single-dish radio data and SMA 
data. Using / = 1.7 x 10 -7 and r c = 2.4 x 10 16 cm in ad- 
dition to the compact component gives a good overall fit 
to the SMA data (Model 5; solid line), except that the 
model provides to much self- absorption at the extreme 
blue-shifted emission. The models are summarized in 
Tabled The two-component model provides the best fit 
to the observed size-velocity relation and the flux den- 
sity profile as shown in Fig. El Fig. El shows that the 
two-component model also fits the available single-dish 



spectra well, with a reduced x 2 of 0.8 (only the total 
integrated intensity of the line was used for this esti- 
mate). Intensities for three additional lines other than 
the ones shown in Fig. El were used in the analysis (see 
iSchoier et alJ l200"fil) . Also, the fit to the central pixel 
spectrum in the CLEANed image is excellent, as shown 
in Fig. 

5. DISCUSSION 

LTE stellar atmosphere models predict that the SiO 
fractional abundance in carbon stars is relatively low , 
typically ~ 5 x 10 -8 (see reviews by lOlasssroldl 119991: 
Millar] I2003L and references therein) about three or- 
ders of magnitude lower than in M-type AGB stars. 
The high SiO fractional abundances, derived in Sect. El 
of 1.5 x 10~ 6 in the inner wind (w3- 8 stellar radii) 
of IRC+10216 can generally not be explained by LTE 
chemistry. 

Departure from LTE could be caused by the vari- 
able nature of AGB stars that induces shocks prop- 
agating through the photosphere thereby affecting its 
chemistry. Models of shocked carbon-rich stellar at- 
mospheres indicate that the SiO fractional abundance 
can be s ignificantly increased by the passage of periodic 
shocks ll Willacv fc Cherchnefj [19981: iHelling fc Winters! 



111 iCher chncfi 2006) . There is also a strong depen- 
dence on the shock strength and possibly this mechanism 
can explain the observed fractional abu ndances, and 
their large spread, of SiO in carbon stars ijSchoier et alJ 
2006) . The lack of including absorption of molecules onto 
dust grains is a shortcoming of current generation chemi- 
cal models and will be required in a full model describing 
the chemistry of AGB stars. In addition, grain surfaces 
could act as catalysts for chemical reactions. 

The model that best explains all available obser- 
vational constraints contains a compact high abun- 
dance component in addition to a more extended low- 
abundance component as illustrated in Sect. El This 



The distribution of SiO in the CSE around IRC +10216 
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is the same behaviour as observed by ISchoier et alJ 
( 2004) for the low-mass-loss-rate M-type stars R Dor and 
L 2 Pup. Schoier et al. introduced a compact compo- 
nent (r e ~ 1 x 10 15 cm) with a high fractional abundance 
(/o ~4 x 10~ 5 ) to explain the small-scale emission. The 
motivation for this is that SiO is expected to effectively 
freeze out onto dust grains in the inner envelope. For a 
high-mass-loss-rate object such as IRC +10216 the freeze 
out is expected t o be very efficient and simple conden- 
sation theory (see lGonzalez Delgado et al.l l)2003f) for de- 
tails) suggests that the SiO fractional abundance should 
be significantly affected within r c w 5 — 10 x 10 14 cm. 

We also find evidence for a sharp decline of the SiO 
fractional abundance close to the photosphere (< 3 stellar 
radii) from a value of w 1.5 x 10 -6 down to «3 x 10~ 8 . 
This lower value is consistent with predictions from 
LTE stellar atmosphere models and could help explain 
why there are no detections of SiO maser emission in 
IRC +10216, or any o ther carbon star, as opposed to Mi- 
type AGB stars (e.g.. lLemne et alJll97^ . 

6. CONCLUSIONS 

Interferometric sub-millimeter SiO J — 5 — > 4 line ob- 
servations performed at the SMA are presented for the 
carbon star IRC+10216. Based on a detailed excita- 
tion analysis and additional IR observations, a fractional 
abundance of SiO in the inner wind of IRC+10216 of 
more than an order of magnitude higher than predicted 
by thermal equilibrium chemistry is derived. A possible 



explanation for the high SiO fractional abundance found 
is a shock-induced chemistry. However, the influence of 
dust grains, both as a source for depletion as well as pro- 
duction of SiO, needs to be further investigated. The 
infrared ro-vibrational data also provide information of 
the wind dynamics close to the star. The velocity of 
the wind has reached within 10% of its terminal value 
already at a distance of 5 x 10 14 cm. 

In addition to the compact high abundance compo- 
nent a spatially more extended region with low SiO abun- 
dance is needed to explain the observations. We find it 
likely that the observed abundance distribution around 
IRC+10216 is the result of SiO depletion in the wind 
at distances > 8 stellar radii. Similar abundance gradi- 
ents have been repo rted for two M-type AGB stars by 
ISchoier et alJ (2004). High spatial resolution interfero- 
metric observations for a larger sample of sources could 
help to solidify this claim. 

Evidence of departure from a smooth wind is found in 
the observed visibilities indicative of density variations a 
factor 2 to 5 on time scales of about 200 years, consistent 
with the densit y modulations found from scatt ered light 
measurements {Mauron fc Hugginslll999l l2000|) . 
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